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Abstract 

Diffractive p meson production has been identified as one of the important processes 
where saturation can be probed at a future Electron-Ion Collider (EIC). A source of un- 
certainty in making predictions for this process lies within the assumed form of the me- 
son light-cone wavefunction. We report here the results of reference [1] where a Regge- 
inspired dipole model was used to extract this wavefunction as well as the corresponding 
leading twist-2 Distribution Amplitude from the current accurate HERA data. In addi- 
tion, we shall check the robustness of the conclusions in reference [1] by using alternative 
Colour Glass Condensate dipole models. 

1. Introduction 

In the dipole model [2, 3], the imaginary part of the amplitude for diffractive p production is written 
as [4] 

QmA x (s, t; Q 2 ) = Y.I d2rd ^V> Q 2 )Kl^ zTe~^N(x, r, A) (1) 
h,h 

where t = — | A| 2 . In a standard notation [1, 4, 5], Vl/?' ^ A and are the light-cone wavefunctions 
of the photon and the p meson respectively while J\f{x, r, A) is the imaginary part of the dipole- 
proton elastic scattering amplitude. The energy dependence of the latter is via the dimensionless 
variable x taken here to be 

x = (Q 2 + 4m})/(Q 2 + s) (2) 

where rrif is a phenomenological light quark mass. 1 Setting t = in equation (1), we obtain the 
forward amplitude used in reference [1], i.e 

9m^ A ( S ,t;Q 2 )| t=0 = S ^ f d\dz^(r,z-,Q 2 )a(x,r)^ h (r,zy (3) 

h,h 

where we have used the optical theorem to introduce the dipole cross-section 

AA(x,r,0) 

a(x, r) = . (4) 
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'Here we shall take m/ =0.14 GeV i.e the value used when extracting the dipole cross-section from F2 data. 



Note that since the momentum transfer A is Fourier conjugate to the impact parameter b, the dipole 
cross-section at a given energy is simply the fe-integrated dipole-proton scattering amplitude, i.e 



a( x , r ) = \ j d 2 bAA(x, r, b) . (5) 
This dipole cross-section can be extracted from the F 2 data since 

F 2 (x,Q 2 )oc J d 2 rdz\V r (r,z;Q 2 )\ 2 a(x,r) (6) 

and the photon's light-cone wavefunctions are known in QED, at least for large Q 2 . The F 2 - 
constrained dipole cross-section can then be used to predict the imaginary part of the forward am- 
plitude for diffractive p production and thus the forward differential cross-section, 



dcr\ 
dt 



= -^(Wl A ( s ,0)) 2 (l + /32), (?) 
t=o 167r 



where /3\ is the ratio of real to imaginary parts of the amplitude and is computed as in reference [1]. 
The t-dependence can be assumed to be the exponential dependence as suggested by experiment [6] : 



dax = da\ 
dt dt 



x exp(-fi|t|) (8) 

=0 



where 



with N = 0.55 GeV 2 . After integrating over t, we can compute the total cross-section a = 
(Tl + eery which is measured at HERA. 2 

Presently several dipole models [7, 8, 9, 10, 11] are able to fit the current HERA F 2 data 
and there is evidence that the data prefer those incorporating some form of saturation [12]. We 
can use the ^-constrained dipole cross-section in order to extract the p light-cone wavefunction 
using the current precise HERA data [6, 13]. This has recently been done in reference [1] using 
the Regge-inspired FSSat dipole model [7] and we shall report the results of this work here. In 
addition, we repeat the analysis using two alternative models [9, 8, 11] both based on the original 
Colour Glass Condensate (CGC) model [14]. They differ from the original CGC model by including 
the contribution of charm quarks when fitting to the F 2 data. Furthermore in one of them [9, 8], the 
anomalous dimension 7 S is treated as an additional free parameter instead of being fixed to its 
LO BFKL value of 0.63. We shall refer to these models as CGC[0.74] and CGC[0.63] models 
where the number in the square brackets stands for the fitted and fixed value of the anomalous 
dimension respectively. For both models, we use the set of fitted parameters given in reference [8]. 
All three models, i.e FSSat, CGC[0.63] and CGC[0.74] account for saturation although in a b- 
(or equivalently t-) independent way. Indeed, at a given energy, the dipole cross-section is equal 
to the forward dipole-proton amplitude given by equation (4) or to the ^-integrated dipole proton 
amplitude given by equation (5). Finally we note that all three dipole models we consider here give 
a good description of the diffractive structure function, i.e F 2 °^ data [15, 16]. 



2 To compare with the HERA data, we take e = 0.98. 



2. Fitting the HERA data 

Previous work [5, 4, 8] has shown that a reasonable assumption for the scalar part of the light-cone 
wavefunction for the p is of the form 

4, B x °(r,z) = A/a 4[z(1 - z)]>»fi^ exp M^A exp (- ^'^A (10) 

^(-^-^) 

and is referred to as the 'Boosted Gaussian' (BG). This wavefunction is a simplified version of 
that proposed originally by Nemchik, Nikolaev, Predazzi and Zakharov [17]. In the original BG 
wavefunction, b\ = 1 while the parameters R\ and M\ are fixed by the leptonic decay width 
constraint and the wavefunction normalization conditions [1]. 

However, when the BG wavefunction is used in conjunction with either the FSSat model or 
any of the CGC models, none of them is able to give a good quantitative agreement with the current 
HERA p-production data. This is illustrated by the large \ 2 values in table 1. As shown in table 2, 
this situation is considerably improved by fitting R\ and b\ to the leptonic decay width and HERA 
data. 3 



Boosted Gaussian predictions 



Dipole model 


X 2 /data point 


FSSat 


310/75 


CGC[0.74] 


262/75 


CGC [0.63] 


401/75 



Table 1: Predictions of the x 2 /data point using the BG wavefunction. 



BG fits 



Model 


X 2 /d.o.f 


FSSat [1] 


82/72 


CGC [0.74] 


64/72 


CGC [0.63] 


83/72 



Table 2: /A.o.f obtained when fitting R\ and b\ to the leptonic decay width and HERA data. 

For the FSSat and CGC[0.63] models, we can further improve the quality of fit by allowing 
for additional end-point enhancement in the transverse wavefunction, i.e. using a scalar wavefunc- 
tion of the form 

Mr, z) = <^ G (r, z) x [1 + c T e + d T t] dD 
where £ = 2z — 1. The results are shown in table 3. 

3 We fit to the same data set and with the same cuts as in reference [1]. 



Improved fits 



Model 


X 2 /d.o.f 


FSSat [1] 


68/70 


CGC[0.63] 


67/70 



Table 3: xVd-o.f obtained when fitting 6a, R\ ct, dr the leptonic decay width and HERA data. 



Best fit parameters 





Rl 




b L 






cLt 


FSSat [1] 


26.76 


27.52 


0.5665 


0.7468 


0.3317 


1.310 


CGC[0.63] 


27.31 


31.92 


0.5522 


0.7289 


1.6927 


2.1457 


CGC[0.74] 


26.67 


21.30 


0.5697 


0.7929 









Table 4: Best fit parameters for each dipole model. 



The best fits obtained with each dipole model are compared to the HERA data in figure 
1, 2 and 3. The corresponding fitted parameters are given in table 4. Note that we achieve a 
lower x 2 / d -°- f = 0.89 with CGC[0.74] than with CGC[0.63] and FSSat for which we obtain 
X 2 /d.o.f = 0.96 and x 2 / d -°- f = °- 97 respectively. Compared to the FSSat and CGC[0.63] fits, 
note that no additional enhancement in the transverse wavef unction is required in the CGC[0.74] 
fit. Nevertheless the extracted wavef unction still exhibits enhancement compared to the old BG 
wavefunction. The extracted light-cone wavefunctions are shown in figure 4. 



3. Distribution Amplitudes 

The leading twist-2 Distribution Amplitude (DA) is given by [1] 

<p{z, n) ~ (l - e"^ 2 ) 2 ) e" m / /A(2)2 [z(l - z)] h ^ , (12) 

where A(z) 2 = 8[z(l — z)] bL /R 2 L . This leading twist DA is only sensitive to the longitudinal 
wavefunction and, as illustrated in figure 5, we expect little variation in the predictions using the 
different dipole models. To compare with existing theoretical predictions for the DA, we compute 
moments, i.e. 

<f%= / dzCp(z^). (13) 
J o 

where by convention [1] 

[ dz¥>(z,M) = l. < 14) 
J o 

In reference [1], we noted that our DA is very slowly varying with \i for fi > 1 GeV, i.e 
our parameterization neglects the perturbatively known /[/-dependence of the DA. This statement 
remains true if we use the CGC[0.63] or CGC[0.74] instead of the FSSat model. 

Our results are compared with the existing predictions in table 5. The moments obtained with 
our best fit, i.e with the CGC[0.74] model, are very similar to those obtained with FSSat model or 
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Fig. 1: Best fits to the HERA total cross-section data. CGC[0.74]: solid; FSSat: dotted; CGC[0.63]: dashed. 

the CGC[0.63]. In all cases, the results are in very good agreement with expectations based on QCD 
sum rules and the lattice. 

Finally, in figure 5 we compare our DAs with that predicted by Ball and Braun [20], at a scale 
p = I GeV. The agreement is reasonable given that in reference [20], the expansion in Gegenbauer 
polynomials is truncated at low order, which is presumably responsible for the local minimum at 
z = 1/2. Certainly all 4 distributions distributions are broader than the asymptotic prediction 

~ Qzil — z). 

4. Conclusions 

We have used the current HERA data on diffractive p production to extract information on the p 
light-cone wavefunction. We find that the corresponding leading twist-2 DA is broader than the 
asymptotic shape and agrees very well with the expectations of QCD sum rules and the lattice. We 
also find that the data prefer a transverse wavefunction with end-point enhancement although the 
degree of such an enhancement is model-dependent. 
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Fig. 2: Best fits to the ZEUS total cross-section data. CGC[0.74]: solid; FSSat: dotted; CGC[0.63]: dashed. 
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Fig. 3: Best fits to the a L /a T data. The HI data are at W = 75 GeV while the ZEUS data are at W = 90 GeV. 
CGC[0.74]: solid; FSSat: dotted; CGC[0.63]: dashed. 




Fig. 5: The extracted leading twist-2 DAs at (i = 1 GeV (CGC[0.74]: solid; FSSat: dotted; CGC[0.63]: dashed) 
compared to the DA of reference [20] also at 1 GeV (long-dashed) and the asymptotic DA (dot-dashed). 



Moments of the leading twist DA at the scale fi 



Reference 


Approach 


Scale /j, 












(This paper) 


CGC[0.74] fit 


~ lGeV 


0.227 


0.105 


0.062 


0.041 


0.029 


(This paper) 


CGC[0.63] fit 


~ lGeV 


0.229 


0.107 


0.063 


0.042 


0.030 


[1] 


FSSat fit 


~ lGeV 


0.227 


0.105 


0.062 


0.041 


0.029 


(This paper) 


Old BG prediction 


~ lGeV 


0.181 


0.071 


0.036 


0.021 


0.014 


[18] 


GenSR 


lGeV 


0.227(7) 


0.095(5) 


0.051(4) 


0.030(2) 


0.020(5) 


[19] 


SR 


lGeV 


0.26 


0.15 








[20] 


SR 


lGeV 


0.26(4) 










[21] 


SR 


lGeV 


0.254 










[22] 


SR 


lGeV 




u - ii=t 0.02 








[23] 


Lattice 


2GeV 


0.24(4) 












6*(1 - z) 


00 


0.2 


0.086 


0.048 


0.030 


0.021 



Table 5: Our extracted values for (£ n ) M , compared to predictions based on the QCD sum rules (SR), Generalised QCD 
Sum Rules (GenSR) or lattice QCD. 
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